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Abstract. We consider multidimensional systems of PDEs of generalized evo- 
lution form with i-derivatives of arbitrary order on the left-hand side and with 
the right-hand side dependent on lower order i-derivatives and arbitrary space 
derivatives. For such systems we find an explicit necessary condition for exis- 
tence of higher conservation laws in terms of the system's symbol. For systems 
that violate this condition we give an effective upper bound on the order of 
conservation laws. Using this result, we completely describe conservation laws 
for viscous transonic equations, for the Brusselator model, and the Belousov- 
Zhabotinskii system. To achieve this, we solve over an arbitrary field the 
matrix equations SA = A*S and SA = —A^S for a quadratic matrix A and 
its transpose A*, which may be of independent interest. 
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1. Introduction 

It is well known that conservation laws are of fundamental importance for clari- 
fying the structure of PDEs. In particular, a common feature of soliton equations 
is to have conservation laws of arbitrarily high order. Existence of higher order 
conservation laws imposes very strong conditions on a system of PDEs. Explicit 
formulation of these conditions would help to classify integrable systems of a given 
type. 

The straightforward study of the conserved current condition is hampered by 
the fact that one is interested in equivalence classes of conserved currents modulo 
trivial ones. Therefore, it is convenient to switch from a conserved current to its 
characteristic, which is the same for equivalent currents and satisfies the equation 
adjoint to the linearization of the initial system [|l|, ^, ^, [l^ . 

Thus a part of the problem is to determine conditions for the adjoint linearized 
equation to have higher order solutions x- In the present article we perform the first 
natural step in this direction. For determined, possibly multidimensional, systems 
of PDEs we find the conditions imposed on the symbol of the system by the fact that 
some higher order vector-functions satisfy the adjoint linearized equation modulo 
lower order terms. 
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More precisely, we consider systems of generalized evolution form 

J, a, 6 = 1, . . . , TO, j = 1, . . . , n, s — 0, . . . ,h — 1, ii + ■ ■ ■ + in ^ N, 

with t-derivatives of fixed order /i > on the left-hand side and with the right-hand 
side dependent on lower order ^-derivatives and arbitrary space derivatives. 

To any m-component vector- function x of the variables Xj, t, and their deriva- 
tives we associate its symbol with respect to the space variables Xj, which is an 
m X m matrix, whose entries are homogeneous polynomials in n variables of degree 
equal to the order o(x) of x with respect to xj. 

Let A be the symbol of the right-hand side of ([^) and S be the symbol of the 
characteristic x of a conservation law for (Q). It turns out that if o{x) > O, where 
O < TV is some constant associated to (^, then the adjoint linearized equation 
implies the matrix equation 

SA = (-1)^+M*S'. (2) 

Here and below ^* is the transpose of A. 

A linear algebra problem arises naturally: for what matrices A does there exist 
a nonzero matrix S such that @) holds? In addition, since for known integrable 
systems there are normally higher conservation laws with nonsingular S, one is 
also interested for which A the matrix S can be taken nonsingular, i.e., when the 
matrices A and ±A* are similar (conjugate). 

In solving these problems there is a difference between the cases n = 1 and n > 1. 
If n = 1, one can switch from homogeneous polynomial in one variable matrices A 
and S to the corresponding matrices of coefhcients and, allowing the coefficients 
to be complex, make use of the Jordan normal form [||. While if n > 1 then the 
entries of the matrices belong to the field of rational functions in several variables, 
which is essentially not algebraicly closed, hence the Jordan normal form is not 
generally applicable. Using more sophisticated algebraic technique, we prove the 
following effective criteria. 

Theorem. For any mx m matrix A with the entries from an arbitrary field F and 
the characteristic polynomial d{X) G F[X] we have the following. 

(1) The matrices A and A* are always similar. 

(2) A nonzero m x m matrix S such that SA ~ —J^S exists if and only if the 
polynomials d{X) and d{—X) have a common divisor of positive degree. 

(3) The matrices A and —A* are similar if and only if all the invariant fac- 
tors di{X) of A (certain divisors of the characteristic polynomial Q) satisfy 
d^{-X) = {-iy°sdd>^)^^(^x). In particular, d{X) = (-l)'^'=e'*(^)d(A), which 
in the case chari^ ^ 2 implies tr A = and, if m is odd, dct A ~ 0. 

Statements 2 and 3 of the theorem give a necessary condition for existence of 
higher conservation laws for systems with odd N + h. In particular, a scalar 
equation (to = 1) of the form (|l|) with odd N + h can not have conservation laws 
of order greater than O. For different ways to write system (|^) in the generalized 
evolution form the symbols A and the resulting conditions are generally different. In 
order to have higher conservation laws a system of PDEs must satisfy all conditions 
obtained from various ways to write it in the generalized evolution form. 
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Let US discuss the previous research on this theme. It seems that only evolu- 
tion systems {h = 1) in one space variable (n = 1) were studied in this respect. 
For such systems equation was obtained by a similar technique in |^ and re- 
discovered in 1^. In j|] it is noticed that SA = —A^S for nonsingular S implies 
det(^ + AI) = ( — l)™det(y4 — AI), which is a weaker version of our Statement 3. 
Here and below I is the unity matrix. In Statement 2 is proved for complex 
matrices, and the corresponding necessary condition for existence of conservation 
laws of order greater than the order of the evolution system is formulated. Even in 
this simplest case our result is stronger, since the upper bound O is normally much 
smaller than the order of the system (see the examples in Section ^ . 

The paper is organized as follows. In Section || the method of characteristics 
of conservation laws is recalled. We specify the method for systems of generalized 
evolution form in Section ^ and derive equation in Section ^. The above theorem 
on quadratic matrices is proved in Section ^. In Section ^ we explicitly formulate 
the obtained necessary conditions for existence of higher conservation laws. Finally, 
Section Q contains some mathematical physics equations of the form (Q), which 
violate these conditions and, therefore, do not have conservation laws of order 
greater than O. This result allows us to describe all conservation laws for two 
basic equations in the theory of viscous transonic gas flows (see, for example, ^ 
^ |ll[ and references therein) and for two popular reaction-diffusion systems: the 
Brusselator model and the Belousov-Zhabotinskii system Section 15.4]. 

2. Characteristics of conservation laws 

This is a brief review of the method of characteristic for computation of conser- 
vation laws. We refer to ^, ^, |l^ for more details. 
Consider a system £ of differential equations 

F,(2/„«^...,i;f,...) = 0, s = l,...,p, (3) 

with independent variables yi, . . . , ya, unknown functions i?-"^, . . . , and 

v) ^ — ■■ -, / = (ii, . . . ,ia) e ^1, 

oVi ■ ■ ■ dya" 

being their derivatives. Here and below 2+ is the set of nonnegative integers and 

\I\=il + --- + la- 

Let J- be the algebra of smooth functions of the variables yi, , and Uj. Although 
the whole set of the variables is infinite, each function is supposed to depend only 
on a finite subset. Denote by J-£ the quotient algebra with respect to the ideal X 
generated by the left-hand sides of equations (H) and their differential consequences 
Dy,^...Dy^^{F,) eJ^. Here 

is 

is the total derivative with respect to yi, where Ij is the multi-index with 1 at 
the i-th place, the other indices of 1^ being zero. For two equivalent functions 
/i, /2 e T, /i - /2 e I, one has fi{yi, v^{yi)) = f2{yi, v'-> {yi)) for any local solution 

v'{y^) to d). 

By definition, the ideal I is invariant under the action of Dy^ , which, there- 
fore, defines a derivation Dy^ of Ts. A conserved current for (||) is an a-tuple 
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J — (Ji, . . . , Ja), where Jk G ■?^£, that satisfies the equation 

a 

Y.DyX'h)=^. (4) 

i=l 

A conserved current is cahed trivial, if it has the form 

jfc = ^A(/:/fc)-I]A(/:fci) 

Kk k<l 

for some functions Cki G ^£ , 1 < fc < ^ < a. Two conserved currents are said to be 
equivalent if they differ by a trivial one. Conservation laws are defined to be the 
equivalent classes of conserved currents. 

Let Jk & J- he such that Jk = Jk + ^- Identity (jj) means that 

a 

Y,DyM) = T.9!Dl{F.) (5) 

i=l sj 

for some functions gj G only a finite number of which is nonzero. Here and 
in what follows for each multi-index / = (zi, . . . ,ia) we denote = Dl\ ■ ■ ■ Di^. 
Consider the functions 

Xs = Y.{-lf^Dl{g\), . = p. (6) 

Generally speaking, representation and functions Xs are not uniquely defined 
by the conserved current J . Assume that system (^) is non-overdetermined and 
nondegenerated |^ , then the corresponding elements Xs = Xs+I J-£ are well- 
defined by J and are all zero if and only if J is trivial. The p-tuple x = {xit ■ ■ tXp) 
is the same for equivalent conserved currents and is called the characteristic (or 
generating Junction 12 ) of the corresponding conservation law. 
In addition, x satisfies the adjoint linearized equation 

K{x) = 0, (7) 
where K is the 6 x p matrix differential operator with the entries 

The homological interpretation of these concepts can be found in |0, ^, |l^ . 

3. Formulas for systems of generalized evolution form 

Consider a system £ of m partial differential equations in n -I- 1 independent 
variables xi, . . . , a;„ and m unknown functions , . . . , u"* of the form 

ah I 

-Q^=F\ ^=l,...,m, (9) 
where are smooth functions of the variables xj, t, u° and the following derivatives 

■ ii, . . . , i„ G Z+, — - — , s < h. (10) 

For each multi-index I — (ii, . . . , i„) e and integer s > we denote 
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Let us describe the algebra Tg. for this system. By induction on s, it follows from 
(H) that the derivatives u\ j with s > h are expressed in terms of 

t, Xj, u\ J, s < h. (12) 

Therefore, for each function of t, Xj, and arbitrary derivatives ( ITTI) there is a 
unique equivalent modulo (^ function of variables ( [l2|) . Thus we can identify J-£ 
with the algebra of smooth functions of variables (|lj). Below all functions are 
supposed to be from JTg. 

The restrictions D^^ , Dt : J-£ of the total derivatives are written in coor- 

dinates (|l^ ) as follows 



k.I.s "^sj 



k,I,s 



A = E -U.^+Y.Diip')o^- (13) 

k,I,s<h-l "^sj k,I "^h-l.J 

The equation DiJq + Y^^=i ^xiJi = for a conserved current J — (Jo, Ji, . . . , J„) 
implies the identity 

A Jo + ^ D.. J, = ^ ^^i^^(,j^^^ _ F'), 

which is the specification of (H) for system According to general formula (||), 
the characteristic x = (Xii • ■ • : Xm) is computed as follows 

X,=Y^i-lf\Di{-^). (14) 

From (^ and (||) we see that the characteristic regarded as a column vector 
satisfies the equation 

{-l)'^D'^{x)=L{x), (15) 
where L is the m x m matrix differential operator with the entries 

[Lh = T.i-^)m o 1^ + eVd^a^ o i^. (16) 

4. Solving the adjoint equation for the highest order terms 

For a (vector-)function / the maximal integer k such that df /du\ / 7^ for some 
Q < s < \ < i < \I\ — k \s called the order of / and denoted by o(/). If 
djjdu\ 7 = for all s, i, / then we set o(/) = —1. The maximal integer s < h 
such that df/dulj ^ for some 1 < i < to, |/| = o(/) is denoted by t(/). The 
order of the characteristic of a conservation law for (^) is called the order of the 
conservation law. 

Consider the ring J-[qi , . . . , qn] of polynomials in n variables with J- as the ring 
of coefficients. For each multi-index / — {ii,.. . ,i„) denote by the monomial 
ql^ • ■ • 9^" G ^[qi, ■•■,'?«]• For any /c-component vector- function x — (xij • ■ • ; Xfe) 
and two integers a, b > let S°^''(x) be the k x m matrix with the entries 

[S"'''(X)]., - E ^n<ll,■■■,qn]■ 



6 SERGEI IGONIN 

We call the nonzero matrix S°'^'^^'*(^^ (x) the symbol of x and denote it by S;^. 

Let A be the symbol of the right-hand side {F^,...,F™) of (^) and set 
N = o(F\ . . . , F™) > 0. By assumption ® , one has t(i^\ . . . , F"') ^ 0. There- 
fore, by definition, 

lAk ^ Y. If,'. (17) 

Applying the Leibniz rule, differential operators ( p^ can be uniquely rewritten 
in the usual form 

[Lh^ E fyDi + Y,g];>Dt. 

\I\<N s<h 

In particular, from definition ( p^ one has 

/l' = (-l)'^7r^ V/:|/|=iV. (18) 

^""■0,1 

We set 

O = max{-l, o(/f ) - TV, 0(5:^") - N}. (19) 
From definition (|l|) it follows that O < N. 

Theorem 1. Let x be the characteristic of a conservation law for (^. If o{x) > O 
then we have 

^^A={-l)'^+^A'S^. (20) 



Proof. Set a — N + o(x), h — t(x) < h. Equation (15) implies, in particular 



1)'^S"'^(A'(X)) = S'''^(L(x)). (21) 



Set 



kj,s<h-l ^'^^J 

Clearly, for any vector- function V' with o(-!/') > and t(?/') < h — 1 we have 

o(A(V.))=o(^), t(A(V')) = t(^A) + 1, Sa(^) = S^. (22) 

By formula (^) and assumption (p^, for any vector- function x with o(x) > 
the vector-function Z)^(x) does not depend on the coordinates j with |/| > a or 
|/| = a, s > b. Moreover, from the whole expression 

only the part 

A''o(^I?^(F'^-)-^)oA'^-i-^(x) 

contributes to ^''^^{D^{x))- 

Therefore, by the definition of S"''', properties (^), and formula (0), we obtain 

S^^\D'i{x)) - S'^'''(a'' o (^i?^(F'=)^-^)oA'-i-''(x)) = 



kJ 



h-l.I 



fc,|/|=o(x) 
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Now let us compute the right-hand side of (|2^). In the case o(x) > O only the 
part J2\i\=N f? of contributes to S"-^{L{x)), since Dt{x) for s < ft- does 
not depend on j with \I\ = a, p > b and the number a is greater than the order 
of any coefhcient or gl^ . Therefore, taking into account formulas (^8|) and (p7[), 
we obtain 

Combining this with ( [2^ ) and one gets ([20|). □ 

In the next section we study the conditions imposed on A by equation (pO|). 

5. Linear algebra problems 

For a ring R we denote by Mk{R) the ring of fc x fc matrices with entries from 
R. Consider an arbitrary field F and denote by F its algebraic closure, i.e., the 
minimal algebraicly closed extention of F. 

Theorem 2. Let A e Mk{F) and let d{X) ~ det(yl — AI) be the characteristic 
polynomial of A. A nonzero matrix S G Mk{F) such that 

SA = -A*S (24) 

exists if and only if the polynomials d{X) and d{—X) have a common divisor of 
positive degree. Equivalently, there are two roots Ai, A2 G F of d{X) such that 
Ai + A2 = 0. 

Proof. The polynomials d{X) and d{^X) have a common divisor of positive degree 
if and only if they have a common root A in F, i.e., both A and —A are roots of 
d(A). 

If (11) holds then for any similar to A matrix A' = CAC^^, C G Mk{F), one 
has S'A' = -A'^S' with S' = C~^*SC~'^. Let us regard A as a matrix from 
Mk{F) D Mk{F). Then we can assume A to be in the Jordan normal form. For 
such A one can easily show that the linear map 

Mu{F) Mk{F), S^SA + A'S, 

has a nontrivial kernel if and only if there are two eigenvalues Ai , A2 G of yl such 
that Ai + A2 = (see || for = C). 

It remains to prove that if (^4|) holds for some S G Mk (F) then there is nonzero 
S' G Mk{F) such that S'A = -A*S". Consider a (possibly infinite) basis {aj of F 
regarded as a vector space over F. One has S = OiSi, where Si G Mk{F). Since 
A G Mk{F), equation (||) imphes S^A = -A^S,. □ 

Recall some criteria for two matrices A, B ^ Mk{F) to be similar (see, for 
example, Chapter 13]). Consider the ring F[X] of polynomials in one variable. 
A matrix C G AIk{F[X]) is said to be unimodular if det C is nonzero and belongs to 
F. For any matrix A G Mk{F) the matrix A — AI G Mk{F[X]) admits a canonical 
decomposition 

A-XI = CiDC2, Ci, £>, C2 GMfc(F[A]), (25) 
such that Ci, C2 are unimodular, while D is diagonal. Moreover, the polynomi- 
als di ~ [D]ii have the leading coefficient 1, and di+i is divisible by di for each 
i = 1, . . . , k — 1. Then the k polynomials di are defined uniquely by A and are 
called the invariant factors of A. Note that there is a simple procedure to compute 
the invariant factors |^ Chapter 13]. 
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Proposition 1 Chapter 13]). Two matrices A, B e Mk{F) are similar if and 
only if they have the same invariant factors. 

We call a polynomial d(A) = 'E.a^X' e F[\] skew if d{-\) = (-l)'i°s'*(^)(i(A), 
i.e., for all i = dcgd{X) + 1 mod 2 one has = 0. 

Theorem 3. For any A e Mk{F) we have the following. 

(1) The matrices A and A* are similar. 

(2) The matrices A and —A* are similar if and only if each invariant factor 
of A is skew. In this case the characteristic polynomial is also skew. In 
particular, in the case char_F ^ 2 we have trA = and, if k is odd, 
det A = 0. 

Remark 1 . Note that in the case char F = 2 the second statement of this theorem 
as well as Theorem |2| are trivial. 

Proof. Consider canonical decomposition ( p5| ) for A. Taking the transpose, we 
obtain 

A* - AI = C2'DCi\ (26) 
which is a canonical decomposition for A*, since Ci*, C2* are clearly unimodular. 
Therefore, the invariant factors of A* are the same, which, by Proposition |l|, implies 
that A and A* are similar. 

Multiplying ( p6| ) by —1 and substituting — A in place of A, we obtain 

-A* - AI = -C2*(-A)i:>(-A)Ci*(-A). (27) 

Denote C[ = -C2*(-A), C'^ = rCi*(-A), and D' = D{-X)T, where T G Mfc(F) is 
the diagonal matrix with the entries [T]ii = (— l)'i°g[^l". From ( [27| ) we obtain the 
canonical decomposition —A* — AI = C[D'C2 for —A*- . According to Proposition ^ 
A and — A* are similar if and only if D' = D, which says that all the invariant factors 
[D\ii of A are skew. In this case the characteristic polynomial is also skew, since 
from ( psi ) it is evidently equal to the product of the invariant factors multiplied by 

(-l)^ □ 

6. Necessary conditions for existence of higher conservation laws 

According to Theorem |^, a necessary condition for existence of conservation laws 
for (|^) of order greater than O is that there is a nonzero m x m matrix with 
entries from T[qi, ...,(?„] such that (|2^) holds. Let us treat A and S;,, as matrices 
with entries from the field F of rational functions in n variables (71, . . . , g„. Then 
Theorem ^ implies the following. 

Theorem 4. If N + h is odd and system possesses a conservation law of order 
greater than O, then the characteristic polynomial d{X) = det(yl — AI) and the 
polynomial d{—X) have a common divisor of positive degree. Equivalently, there are 
eigenvalues Ai, A2 of A {possibly in the algebraic closure of F) such that A1+A2 = 0. 

Remark 2. Evidently, introducing the new dependent variables 

u'''' = i = 1, . . . , m, s = 0, . . . , ft. - 1, 

we can rewrite (^) in the usual evolution form. But if ft > 1 then the symbol of 
the right-hand side of the obtained evolution system has zero determinant and, 
therefore, automatically satisfies the condition in Theorem ^, even if the initial 
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system does not meet this condition. Therefore, it is essential to consider the 
generahzed evolution form. 

Analyzing examples of known soliton equations, we can conjecture that for (j^) 
to be integrable there must exist higher order conservation laws with nonsingular 
matrix S^. Therefore, it is worth formulating a necessary condition for existence 
of such nonsingular conservation laws. According to Theorem |^ (|^), we obtain the 
following criterion. 

Theorem 5. If N + h is odd and system ^ has a nonsingular conservation law 
of order greater than O then all the invariant factors of the matrix A are skew. In 
this case its characteristic polynomial is also skew. In particular, tv A ^ and, if 
m is odd, det A — 0. 

If TO — 1 then A is a nonzero 1x1 matrix, and Theorem |^ implies the following. 

Theorem 6. A scalar equation (to — 1) of the form (|^) with odd N + h can not 
have conservation laws of order greater than O. 

Remark 3. To obtain stronger conditions, it is sometimes useful to write a system of 
PDEs in several ways in the form (||). For example, for a scalar evolution equation 
Ut = Uxxx + Uyy + f(u, Ux, Uy , Uxx) thc couditiou is empty, since the sum of h = 1 
and A'^ = 3 is even. But rewriting the equation in the generalized evolution form 
with respect to y as follows Uyy = —Uxxx + ut ~ f{u, Ux, Uy, Uxx), we see, according 
to Theorem ^, that there are no higher conservation laws. 

Remark 4. If A^+/i is even then, according to Theorem || (1), equation ( |2^ ) is always 
solvable and imposes no restrictions on the symbol A of (£). That is, equation (15) 
is solvable with respect to the highest order terms. In this case, to obtain nontrivial 
conditions for (H) to have higher conservation laws, deeper analysis of (^5|) is needed. 

7. Examples 

7.1. The viscous transonic equation. The nonlinear viscous transonic equation 

a 

Utt = -Uxxx + UxUxx - —ut (28) 

describes the asymptotic form of a gas flow in the sonic region (see ^| and 
references therein). The following conserved currents for ( p8[ ) were found in |^ 

(utt", - ^t"), {utt+{ci-l)u,uxxt~^t). (29) 

All other conserved currents mentioned in Q are trivial. 

Let us show that (^9|) span the whole space of conservation laws for (psj). We 
have 

- = Dl + UxDl + UxxDx + -Dt - ^. 
According to (|l^, one has O = — 1. By Theorem |[ since N + h = 5 is odd, the 



L = Dl + Dloux - DxOUxx + Dto^ = Dl + UxDI + UxxDx + t-D* 



characteristic x of conservation law is a function of x, t only. Equation (15) 
reads 

a a 

Xtt — Xxxx ~t" UxXxx ~t" llxxXx ~t" ~i^Xt 



This implies x^: = and 



Xtt = jXt - -^X- (30) 
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The characteristics t", t of conserved currents (g9|) span the space of solutions to 
(^0|). Therefore, ( [2S| ) span the space of conservation laws. 

7.2. Another equation for viscous transonic flows. The equation 

Vyy = 2Vxt + VxVxX - ^ZZ - fJ'VxXX, (3l) 

where is a nonzero real constant, models nonstationary transonic flows around a 
thin body with effects of viscosity and heat conductivity when the velocity of the 
gas is close to the local speed of sound, see ^ |ll| and references therein. 
For this equation we obtain 

By definition (|l9|), we have O = — 1. According to Theorem ^, since N + h = b 
is odd, the characteristic x of ^J^Y conservation law is a function of x, y, z, t only. 
Equation (^5|) implies 

Xx 0, Xyy + Xzz = 0. (32) 



Each function x satisfying (32) is indeed the characteristic of the conserved current 

Dy{xvy - Xyv) + D^{xvz - Xzv) + D^{iMxVxx - ^xvl - 2xvt) = 0. (33) 

Therefore, any conserved current for ( ^l|) is equivalent to a conserved current of 
the form (^). 

7.3. The Brusselator model. The Brusselator model governing certain chemical 
reactions is the following multidimensional system 

Vt = ^-^ + v'^w - (b + l)v + a, 

(34) 

——2 V W + OV, 

z— 1 

where a, 6, and c ^ are real parameters [|l0[ Section 15.4]. By definition ( |l6|) we 
have 

J2i + 2vw - (b + 1) -2v_w + b 



L 



The symbol A is diagonal with [A]ii — J2i^ ^^"^ [A]22 — cJ2i1i- By defini- 
tion i\m, we have O = —1. By Theorem ^ since [A]ii ^ 0, [A]22 7^ 0, and 
+ [A]22 ^ 0, the characteristic x = (xijXz) of any conservation law is a 
function of Xi, t only. For such x equation (|l5| ) reads 

= E ^ + - {b + l))xi - {2vw - b)x2, 

(35) 



dX2 2 , 9^X2 2 
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Evidently, ( pS]) implies xi — X'2- Then (|35|) becomes 

Xi, 



dt ^ dxi 

' (36) 



9X1 ^ 9 Xl 

' dt dx^ ' 



Clearly, for c = 1 the only solution to ( |36D is xi = Oi while for c 7^ 1 the general 
solution is 

ct 

Xi = X2 = G(a;,) exp( -), (37) 

c — 1 

where G{xi) is an arbitrary solution to the equation 

G+(^-l)E^ = 0- (38) 

i * 

Each solution ( |37| ) is indeed the characteristic of the conserved current 

A(Gexp(^)(« + u;)) + 
c — 1 



E 

i=\ 



Thus if c 7^ 1 then these conserved currents span the space of conservation laws; 
while if c = 1 then there are no nontrivial conservation laws for (jsj). 

7.4. The Belousov-Zhabotinskii system. This system describes certain chem- 
ical reactions and reads Section 15.4] 

——2 + v{\ — V — rw) + Lrw, 

I '■ 

-—^ — bvw — Mw. 

iJX: 

Here r, L, 6, M are real constants. Evidently, Theorem ^ implies that this system 
does not possess conservation laws of nonnegative order. Similarly to the above 
examples, analysis of equation (p^ ) for characteristics of order —1 shows that in 
the nonlinear case 6^0 there are no nontrivial conservation laws at all. 
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